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ABSTRACT 

We  present  the  results  of  an  experiment  to  image  the  interacting  binary  star  ft  Lyrae  with  data  from  the  Navy 
Prototype  Optical  Interferometer  using  a  differential  phase  technique  to  correct  for  the  effects  of  the  instrument 
and  atmosphere  on  the  interferometer  phases.  We  take  advantage  of  the  fact  that  the  visual  primary  of  /J  Lyrae  and 
the  visibility  calibrator  we  used  are  both  nearly  unresolved  and  nearly  centrally  symmetric,  and  consequently  have 
interferometric  phases  near  zero.  We  used  this  property  to  correct  for  the  effects  of  the  instrument  and  atmosphere  on 
the  phases  of  /3  Lyrae  and  to  obtain  differential  phases  in  the  channel  containing  the  Her  emission  line.  Combining 
the  Hcr-channel  phases  with  information  about  the  line  strength,  we  recovered  complex  visibilities  and  imaged  the 
Ha  emission  using  standard  radio  interferometry  methods.  Our  images  show  the  position  of  the  Ha-emitting  regions 
relative  to  the  continuum  photocenter  as  a  function  of  orbital  phase,  indicating  a  major  axis  line  of  nodes  along 
il  —  249°  ±4°.  The  orbit  is  smaller  than  previously  predicted,  a  discrepancy  that  can  be  alleviated  if  we  assume  that 
the  system  is  at  a  larger  distance,  or  if  the  stellar  continuum  contribution  to  the  Ha  channel  was  underestimated.  We 
do  not  detect  a  jet  in  the  Ha  images,  which  may  be  due  to  the  limited  resolution  of  the  observations  along  the  direction 
perpendicular  to  the  orbital  plane.  We  find  that  the  differential  phase  results  are  consistent  with  those  obtained  from 
a  more  standard  analysis  using  squared  visibilities  (V2s)  and  closure  phases,  which  also  indicate  an  Ha  disk  radius 
of  0.6  ±0.1  mas,  and  AT  =  1.30  ±0.1  and  A R  =  1.20  ±0.1  mag  for  the  magnitude  difference  between  the  stars. 

Key  words:  binaries:  eclipsing  -  methods:  data  analysis  -  stars:  imaging  -  stars:  individual  (ft  Lyrae)  -  techniques: 
high  angular  resolution  -  techniques:  interferometric 


1.  INTRODUCTION 

One  of  the  major  challenges  in  optical  interferometry  is 
recovering  images,  since  fringe  visibility  phases  are  badly 
corrupted  by  the  atmosphere  on  short  timescales.  A  few  attempts 
have  been  made  to  reconstruct  images  from  the  squared  visibility 
( V 2),  an  unbiased  estimator  of  the  visibility  amplitude  (e.g., 
Baldwin  et  al.  1996;  Quirrenbach  et  al.  1994),  but  in  most  cases, 
images  in  optical  interferometry  are  produced  from  techniques 
that  use  V2  and  closure  phases  (e.g.,  Buscher  1994;  Ireland  et  al. 
2006;  Monnier  et  al.  2007;  Thiebaut  2008).  These  techniques 
have  been  tested  and  shown  to  reliably  reconstruct  images 
(Lawson,  et  al.  2006;  Cotton  et  al.  2008).  However,  since  a 
closure  phase  is  the  sum  of  visibility  phases  around  a  triangle  of 
interferometer  array  elements,  and  the  atmospheric  turbulence 
effects  cancel  in  the  sum,  only  a  fraction  1  —  2/N  of  the 
phase  information  from  an  array  of  N  elements  is  recovered. 
Therefore,  a  challenge  for  optical  interferometric  imaging  is 
to  recover  as  much  of  the  remaining  phase  information  as 
possible. 

Here,  we  describe  the  technique  we  developed  to  obtain  Ha 
differential  phases  from  Navy  Prototype  Optical  Interferometer 
(NPOI)  observations  of  the  interacting  binary  star  /3  Lyrae. 
We  were  able  to  recover  all  baseline  phase  information  for  the 
Ha  channel  of  this  source.  This  allowed  us  to  obtain  complex 
visibilities,  and,  for  the  first  time,  to  obtain  Ha  images  using 
standard  radio  interferometric  reconstruction  techniques. 


Differential  phases,  sometimes  referred  to  as  phase  referenc¬ 
ing  (see  Quirrenbach  2000,  Monnier  2003  for  a  discussion  on 
the  subject)  is  an  interesting  and  powerful  technique.  Here  one 
needs  multiwavelength  observations  obtained  simultaneously 
and  a  priori  knowledge  about  the  structure  of  the  source  in  a 
region  of  the  observed  spectrum.  From  that  knowledge,  one 
calculates  the  expected  phases  in  that  spectral  region,  corrects 
for  the  atmospheric  and  instrumental  effects,  and  interpolates 
or  extrapolates  the  corrections  to  the  wavelengths  of  interest. 
One  of  the  first  attempts  to  use  this  technique  was  presented 
by  Vakili  et  al.  (1997,  1998),  who  used  Ha  and  He  i  667.8 
nm  GUT  observations  of  P  Cygni  and  £  Tauri  to  detect  struc¬ 
tures  in  the  line-emitting  regions  of  these  stars.  However,  they 
did  not  reconstruct  images  based  on  their  measurements;  their 
conclusions  were  based  only  on  the  analysis  of  the  phases  and 
visibilities  of  the  stars.  Another  application  of  the  differential 
phase  technique  was  presented  by  Akeson  et  al.  (2000),  who 
used  the  Palomar  Testbed  Interferometer  to  detect  variations 
of  phase  versus  wavelength  in  data  from  the  binary  star  i  Peg. 
An  ideal  class  of  objects  to  which  the  differential  phase  tech¬ 
nique  can  be  applied  is  Be  stars.  These  sources  are  composed 
of  a  B  star,  where  the  photosphere  radius  /?*  <  10  RQ  (cor¬ 
responding  to  <  1  mas  for  nearby  sources)  usually  is  unre¬ 
solved  by  interferometric  observations,  surrounded  by  a  disk 
of  gas  with  a  diameter  of  ~  100  RQ  (a  few  mas  for  nearby 
sources),  which  is  detected  in  Ha  line  emission  (Tycner  et  al. 
2005). 
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Table  1 

Observations  Log 


Date 

UT  Range 

Julian 

Date 

No.  of 

Scans 

Baselines 

Hour  Angle 

Orbital 

Phase 

2005  May  09 

9.63.. 

.11.90 

3499.95833 

5 

E6-AC,  AW-E6,  AW-AC,  AE-AC,  AN-AC,  AE-AN 

-1.50. 

.0.79 

0.47 

2005  May  18 

9.85.. 

.11.80 

3508.95833 

10 

E6-AC,  AW-E6,  AW-AC,  AE-AC 

-0.67. 

.1.28 

0.17 

2005  May  19 

9.58.. 

.11.75 

3509.95833 

9 

E6-AC,  AW-E6,  AW-AC,  AE-AC,  AN-AC,  AE-AN 

-0.89. 

.1.29 

0.24 

2005  May  24 

9.08. . 

.11.63 

3514.95833 

8 

E6-AC,  AW-E6,  AW-AC,  AE-AC,  AN-AC,  AE-AN 

-1.10. 

.1.48 

0.63 

2005  May  26 

8.98.. 

.11.36 

3516.95833 

10 

E6-AC,  AW-E6,  AW-AC,  AE-AC,  AN-AC,  AE-AN 

-1.04. 

.1.34 

0.78 

Notes.  Column  3  shows  the  Epoch  of  the  observations  (HJD— 2450000)  calculated  for  UT  =  1 1:00  hours,  roughly  the  middle  of 
the  observations.  The  maximum  lengths  of  the  baselines  shown  in  Column  5  are  E6-AC  34.3  m,  AW-E6  53.3  m,  AW-AC  22.2  m, 
AE-AC  18.9  m,  AN-AC  22.9  m,  and  AE-AN  34.9  m.  We  always  used  AC  as  the  reference  station. 


The  eclipsing  binary  star  /i  Lyrae  (HD  174638,  HR  7106, 
FK5  705)  has  been  extensively  studied  since  it  was  discovered 
to  be  an  eclipsing  system  over  220  years  ago  (Goodricke  & 
Englefield  1785).  Nevertheless,  it  remains  one  of  the  most 
baffling  binary  systems  known  (see  Harmanec  2002  for  a 
comprehensive  review).  The  current  model  consists  of  a  ~  3  MG 
B6-8II  star  (star  2  in  Harmanec’s  notation),  which  has  filled  its 
Roche  lobe  and  is  losing  mass  at  a  rate  of  2  x  10-5  MG  year-1  to 
a  ~  13  Mg  early  B  star  (star  1),  which  is  completely  embedded 
in  and  hidden  by  a  thick  accretion  disk.  This  accretion  disk 
simulates  a  pseudophotosphere,  which  appears  to  an  observer 
to  have  a  spectrum  of  an  A5III  star  (for  a  diagram  of  the  P  Lyrae 
system  we  refer  to  Figure  8  in  Linnell  et  al.  1998  and  Figure  8  in 
Hoffman  et  al.  1998).  The  orbit  is  very  nearly  circular,  the  orbital 
period  is  a  little  over  12.9  days,  and  the  period  is  increasing  by 
19  s  year-1.  The  system  lies  at  a  distance  of  270  ±  39  pc,  as 
measured  by  Hipparcos  (Perryman  et  al.  1997).  The  spectrum 
of  /3  Lyrae  contains  at  least  six  spectroscopic  line  components 
(Harmanec  2002).  The  three  major  components  are  (1)  strong 
absorption  lines  from  star  2  used  to  derive  the  radial  velocity 
curve  having  an  amplitude  of  370  km  s  ,  (2)  “shell”  absorption 
components  seen  in  the  hydrogen,  helium,  and  many  metallic 
lines,  believed  to  arise  from  the  accretion  disk,  and  (3)  strong 
lines  of  hydrogen  and  helium  seen  in  emission,  which  vary  with 
the  orbital  phase  and  time.  One  of  the  features  that  makes  f 
Lyrae  unique  among  contact  binaries  is  the  possibility  that  the 
strong  emission  lines  come  from  a  bipolar  jet  that  has  formed 
near  the  point  where  material  streaming  from  star  2  interacts 
with  the  accretion  disk  around  star  1  (Harmanec  et  al.  1996). 

This  paper  is  organized  in  the  following  way.  In  Section  2, 
we  present  the  observations  and  reductions.  In  Section  3,  we 
describe  the  corrections  applied  to  the  data  in  order  to  obtain 
differential  phases  and  describe  the  results  obtained  for  p  Lyrae. 
In  Section  4,  we  describe  the  imaging  of  the  Ha  emission 
of  this  source.  In  Section  5,  we  obtain  information  about  the 
orbit  of  this  system  from  the  displacement  of  the  Ha  emission 
relative  to  the  continuum  photocenter.  In  Section  6,  we  compare 
the  results  obtained  from  the  analysis  of  the  Ha  images  with 
those  obtained  using  a  traditional  V 2  analysis,  and  in  Section  7, 
we  give  a  summary  of  our  results. 

2.  OBSERVATIONS  AND  REDUCTIONS 
2.1.  Interferometry 

Our  observations  were  taken  with  the  NPOI  (Armstrong  et  al. 
1998)  on  five  different  nights  between  2005  May  9  and  2005 
May  26.  The  observations  were  made  with  two  spectrographs, 
simultaneously  recording  fringes  in  16  spectral  channels  in  the 
wavelength  range  560-860  nm.  Between  one  and  three  baselines 
were  measured  with  each  spectrograph.  The  maximum  baseline 


lengths  ranged  from  18.9  to  53.3  m.  In  Table  1,  we  present  a  log 
of  the  observations,  including  the  number  of  scans  observed  each 
night,  the  baselines  used,  and  the  range  of  hour  angles  covered 
by  the  observations.  Figure  1  presents  the  ( u ,  u)-plane  coverage 
achieved  by  the  observations  on  each  night.  For  reference, 
this  figure  also  shows  the  array  layout  with  the  stations  and 
baselines  used.  The  observations  of  f>  Lyrae  were  interleaved 
with  observations  of  a  nearby  calibrator  star  y  Lyrae  (B9III, 
V  —  3.24  mag),  which  has  a  similar  spectral  type  and  magnitude 
to  that  of  p  Lyrae.  Each  calibrator  or  source  scan  had  a  duration 
of  30  s. 

Throughout  this  paper  we  will  use  orbital  phases  of  P  Lyrae 
based  on  the  latest  ephemeris  of  primary  eclipses  from  Rreiner 
(2004). 9  We  assumed  a  single  orbital  phase  for  each  night, 
corresponding  to  the  midpoint  of  the  observations  (Table  1). 
This  assumption  will  have  no  effect  on  the  final  results  of  this 
paper,  because  our  observations  usually  lasted  for  ~  2  hr,  which 
corresponds  to  ~  0.6%  of  an  orbital  period.  Comparing  the 
values  obtained  from  this  ephemeris  with  those  from  Harmanec 
&  Scholz  (1993),  we  find  that  the  two  predictions  agree  at  the 
2%  level,  which  is  sufficient  for  the  purposes  of  the  analysis 
presented  here.  As  a  final  check,  we  compared  the  predicted 
times  of  minima  with  visual  measurements  from  AAVSO,10 
spread  over  a  period  of  two  months  around  our  observations. 
We  used  only  validated  observations  available  at  the  AAVSO 
site  and  got  an  agreement  better  than  0.25  days  between  the 
observed  and  predicted  minima,  which  corresponds  to  ~  2%  of 
the  orbital  period. 

The  initial  data  reductions  followed  the  incoherent  and  co¬ 
herent  integration  methods  described  in  Peterson  et  al.  (2006a). 
In  the  case  of  the  incoherent  integration,  we  start  with  data  that 
were  processed  to  produce  squared  visibilities  ( V 2)  averaged 
into  1  s  intervals.  These  data  are  then  flagged  to  eliminate  points 
with  fringe  tracking  and  pointing  problems  (Hummel  et  al.  1998, 
2003a).  The  flagged  data  set  is  then  used  to  bias  correct  the  V 2 
data,  which  are  then  calibrated.  The  calibrated  V2  values  are 
later  used  to  estimate  the  amplitudes  of  the  complex  visibilities. 

For  the  coherent  integration  of  complex  visibilities,  we  use 
the  technique  presented  by  Hummel  et  al.  (2003b;  see  Jorgensen 
et  al.  2007  for  a  different  algorithm).  We  treat  the  visibilities  in 
each  2  ms  frame  (the  instrumental  integration  time)  as  phasors. 
To  determine  the  phase  shift  from  one  frame  to  the  next,  we 
compute  the  power  spectrum  of  the  channeled  visibility  as  a 
function  of  delay  and  measure  the  phase  difference  between  the 
peak  of  the  fringe  envelope  and  the  nearest  fringe  peak.  We 
rotate  the  phasors  to  correct  for  this  residual  phase  and  average 
the  complex  visibilities  in  200  ms  subscans,  which  are  used 


9  http://www.as.wsp.krakow.pl/ephem/LYR.HTM 

10  www.aavso.org 
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Figure  1.  u  —  v  coverage  of  the  different  nights.  The  small  points  show  the  continuum  channels,  while  the  large  ones  show  Ha?.  The  observing  date  is  presented  in  the 
top  right  corner  of  each  panel.  The  bottom  right  comer  of  the  figure  shows  the  array  layout.  The  astrometric  stations  are  shown  as  black  squares  and  E06  is  shown  as 
a  black  dot.  We  also  show  the  baselines  used  on  all  nights  as  solid  lines  and  the  baselines  that  were  used  only  on  some  nights,  those  involving  AN,  as  dashed  lines. 
The  maximum  baseline  lengths  can  be  found  in  Table  1 . 


to  determine  the  visibility  phases.  In  Section  3,  we  describe 
other  corrections  applied  to  these  data  in  order  to  obtain  Ha 
differential  phases  for  /3  Lyrae. 

2.2.  Spectroscopy 

To  estimate  the  strength  of  the  Ha  emission  in  the  p  Lyrae 
system,  we  have  obtained  a  high-resolution  spectrum  using  a 
fiber-fed  Echelle  spectrograph  at  the  Lowell  Observatory’s  John 
S.  Hall  telescope.  The  spectrum  was  acquired  on  2005  May 
19,  which  corresponds  to  the  middle  of  our  interferometric 
run.  The  spectroscopic  data  were  processed  using  standard 
routines,  developed  specifically  for  the  instrument  used  in  the 
observations  (Hall  et  al.  1994).  The  final  reduced  spectrum  in 
the  Ha  region  reaches  a  resolving  power  of  10,000  and  a  signal- 
to-noise  ratio  (S/N)  of  few  hundreds  (see  Figure  2). 

We  obtain  an  equivalent  width  (EW)  of  —1.5  ±  0.05  nm  for 
the  Ha  profile  shown  in  Figure  2.  Most  of  the  uncertainty  in 
the  EW  measurement  is  due  to  the  uncertainty  associated  with 
continuum  normalization,  which  is  estimated  to  be  of  the  order 
of  3%.  To  obtain  the  total  emission  in  the  Ha  line,  we  also 
need  to  correct  for  the  underlying  stellar  absorption  line  that  has 
been  filled  in  by  the  emission  from  the  circumstellar  material. 
We  estimate  the  EW  of  the  underlying  absorption  line  using 
the  calibration  developed  by  Cote  &  Waters  (1987).  Using  a 
reddening-corrected  color  ( B  —  V)o  =  —0.05  mag  for  ft  Lyrae 
from  Dobias  &  Plavec  (1985),  we  obtain  an  EW  of  0.72  nm 
for  the  Ha  absorption  component,  and  therefore  we  estimate 
that  the  EW  of  the  total  Ha  emission  is  —2.2  nm.  Considering 
that  the  NPOI  Ha  channel  has  a  width  of  15  ±  1  nm  and  the 
total  emission  line  width  is  —2.2  nm,  we  derive  a  fractional 
contribution  from  the  stellar  photosphere  (ccnt)  to  the  total  flux 
in  the  channel  of  0.87  ±  0.01. 

3.  DIFFERENTIAL  PHASES 

The  challenge  of  recovering  phase  information  from  opti¬ 
cal  interferometric  observations  results  from  atmospheric  and 
instrumental  effects,  changing  the  phase  calibration  on  a 
timescale  of  a  few  milliseconds.  Fringe  tracking  can  reduce 


Figure  2.  High-resolution  Ha?  profile  of  ft  Lyrae,  normalized  to  the  mean 
continuum  level. 


these  variations  but  not  to  the  level  needed  for  good  imaging. 
Even  with  fringe  tracking,  the  solution  starts  with  ( 1 )  short  expo¬ 
sures  that  freeze  the  fringe  motion,  (2)  software  that  measures 
and  removes  these  variations  and,  (3)  a  coherent  integration 
of  hundreds  of  milliseconds  (or  more)  needed  to  increase  the 
S/N  to  a  level  sufficient  for  the  rest  of  the  data  processing. 
The  result  is  a  measurement  of  the  real  and  imaginary  parts 
of  the  visibility.  The  uncorrected  phase  variations  during  the 
coherent  integration  consist  of  a  constant  part  and  zero  mean 
fluctuations.  The  fluctuations  reduce  the  measured  visibility  am¬ 
plitude  but  have  no  effect  on  the  phase.  Thus,  these  visibility  data 
can  be  processed  using  closure  relationships  or  self-calibration. 
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Figure  3.  Coherently  integrated  phases  of  three  scans  of  the  calibrator  y  Lyrae,  observed  on  2005  May  9  with  baseline  AW0-E06  at  UT  =  9:38  (left),  1 1:23  (middle), 
and  1 1 :55  (right).  Each  line  corresponds  to  a  different  200  ms  subscan. 


However,  with  a  small  number  of  stations,  there  are  signifi¬ 
cantly  more  baselines  than  closure  phases,  so  a  significant  frac¬ 
tion  of  phase  information  is  lost.  By  contrast,  the  differential 
phase  technique  described  here  provides  good  phase  informa¬ 
tion  for  every  baseline.  The  idea  is  to  recognize  that  the  observed 
phase  consists  of  three  terms,  viz.,  the  desired  source  phase  0o, 
and  offsets  caused  by  the  atmosphere  and  the  instrument: 

0obs  =  00  +  0atm  "b  0inst-  (1) 

All  terms  have  implicit  wavelength  dependences.  The  differ¬ 
ential  phase  technique  consists  of  (1)  measuring  0obs  at  several 
wavelengths,  (2)  modeling  the  three  right-hand  side  terms  in 
Equation  (1),  and  (3)  deriving  parameters  of  the  model  by  fit¬ 
ting  it  to  the  measured  0O bs.  For  the  NPOI,  0jnst  is  stable  and  can 
be  estimated  by  observing  calibration  stars;  no  parameters  need 
to  be  fitted.  The  atmospheric  term,  0atm,  is  given  by 

0atm  =  a(n  -  l)k  +  elk  +  0C,  (2) 

where  a  is  the  differential  air  path  between  two  beams,  n  is  the 
refractive  index  of  air,  k  the  wavenumber  (2tt /k),  d  the  vacuum 
path  difference,  and  0C  is  a  wavelength-independent  phase  offset 
(Jorgensen  et  al.  2006).  Three  quantities  need  to  be  fitted:  a, 
d,  and  0C.  Thus,  with  at  least  four  spectral  channels,  some 
information  about  the  source  can  be  obtained  for  each  baseline. 

3.1.  Instrumental  and  Atmospheric  Corrections 

The  Ha  emitting  circumstellar  disks  of  Be  stars  constitute  one 
of  the  best  cases  for  which  one  can  calculate  instrumental  and 
atmospheric  corrections  to  the  visibility  phases.  This  correction 
is  based  on  observations  of  calibrator  stars  and  the  continuum 
channels  of  the  star  of  interest.  Since  the  photospheres  of  these 
stars  usually  do  not  deviate  significantly  from  point  symmetry, 
we  can  model  the  continuum  as  having  zero  intrinsic  phase. 
Even  in  the  case  of  rapidly  rotating  stars,  in  which  the  stellar 
image  may  deviate  from  point  symmetry,  the  phase  should 
vary  smoothly  and  slowly  with  wavelength  (Yoon  et  al.  2007; 
Peterson  et  al.  2006a,  2006b).  The  same  principle  applies 
to  calibration  stars.  These  properties  allow  one  to  use  the 
information  from  the  continuum  channels  to  calculate  the 
appropriate  corrections  to  the  phases  and  interpolate  them  over 
the  channel  containing  the  Ha  line. 

A  first  example  of  the  instmmental  and  atmospheric  effects  on 
the  phases  of  a  partially  resolved  star  is  shown  in  Figure  3.  Each 
panel  of  this  figure  shows  a  different  scan  of  y  Lyrae,  while 
the  different  lines  in  each  panel  represent  individual  200  ms 
subscans.  Note  that  the  lines  of  the  different  subscans  follow  a 
similar  cubic  like  shape  as  a  function  of  wavenumber,  with  some 


scatter.  This  cubic  curve  can  be  explained  as  an  instrumental 
phase  due  to  uncompensated  amounts  of  glass  along  the  different 
beams.  This  curve  does  not  change  significantly  from  scan  to 
scan,  or  from  night  to  night.  A  further  inspection  of  this  figure 
shows  that,  besides  this  cubic  component,  one  can  also  see 
a  significant  variation  of  the  phase  among  different  subscans. 
These  variations  are  due  to  the  differential  air  path  between  the 
two  beams,  and  the  shape  of  the  curve  is  due  to  the  fact  that 
(n  —  1)  has  a  significant  dependence  on  k2  (Owens  1967).  One 
can  also  see  in  Figure  3  that  the  scatter  is  much  smaller  in  the  left 
panel,  during  the  first  scan  of  the  night,  increasing  significantly 
toward  later  times  (right),  consistent  with  the  worsening  of  the 
seeing  conditions  through  the  night. 

In  order  to  obtain  differential  phases  of  Be  stars,  one  needs 
to  correct  the  data  for  the  two  components  described  above. 
For  the  instrumental  phases  we  use  the  values  determined  by 
averaging  all  the  200  ms  subscans  of  calibrator  stars  observed 
throughout  the  night.  Since  the  effects  of  the  differential  air  path 
and  vacuum  delay  are  random  and  constitute  an  additive  term  to 
the  intrinsic  phases,  one  can  reasonably  assume  that  they  should 
average  down  to  zero  over  the  course  of  a  night.  Also,  since 
calibration  stars  are  expected  to  have  zero  phase  in  the  spectral 
range  covered  by  our  observations,  by  averaging  the  complex 
visibilities  of  the  calibrators  one  is  left  only  with  instmmental 
phases,  which  can  then  be  subtracted  from  the  program  star 
visibilities. 

We  demonstrate  the  procedure  used  to  obtain  instrumental 
phases  in  Figure  4.  Here  we  show  the  phases  obtained  by 
averaging  the  visibilities  of  all  720  subscans  of  y  Lyrae  observed 
on  the  same  night  (white  line).  We  can  also  show  that  the 
other  component  affecting  the  phases  of  the  subscans  is  due 
to  differential  air  path  between  the  two  beams.  Following  the 
derivation  presented  by  Armstrong  et  al.  (1998),  we  expand 
Equation  (2)  into  a  power  series  around  k{)  and  rewrite  Equation 
(1)  in  the  following  form: 

0obs  =  00  +  0inst  +  <Y  +  C\k  +  C2^“ .  (3) 

Starting  from  this  equation  and  assuming  that  0o  =  0,  one  can 
then  solve  for  the  coefficients  Co,  c i,  and  C2  for  each  subscan, 
such  that  the  residuals  produce  the  best  fit  to  the  atmosphere 
phases  (0obs  —  0inst).  We  show  in  Figure  4,  as  black  lines,  all 
subscans  after  subtracting  only  the  atmosphere  effects  from  the 
observed  phases.  Note  that  these  lines  follow  the  instrumental 
one.  Figure  4  also  shows  histograms  with  the  distribution 
of  phases  after  the  subtraction  of  atmospheric  effects.  These 
histograms  show  that  the  phases  have  a  much  lower  scatter  in 
the  red  part  of  the  spectrum,  where  the  S/N  is  higher.  The 
uncertainties  in  the  mean  values  of  the  instmmental  phases 
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Figure  4.  Top  panel  shows  the  comparison  between  the  phases  of  individual 
subscans  of  y  Lyrae  from  Figure  3,  corrected  for  atmospheric  and  vacuum 
delay  effects  but  not  for  the  instrumental  phase,  which  are  shown  as  a  white 
line.  The  instrumental  phases  were  obtained  by  averaging  all  subscans  of  this 
calibrator  star  observed  in  this  particular  night.  This  figure  demonstrates  that 
the  average  phase  of  y  Lyrae  is  a  good  estimator  of  the  instrumental  phase.  The 
bottom  panels  show  histograms  of  the  distribution  of  the  corrected  phases  in 
three  different  channels.  The  uncertainty  in  the  mean  value  of  the  instrumental 
phase  depends  on  AT-1/2  (number  of  coherently  integrated  subscans),  which 
in  this  case  is  720.  It  grows  from  0°17  in  channel  1  (reddest  one),  to  0°24  in 
channel  8,  and  0°42  in  channel  16  (bluest  one). 


depend  on  N  1/2  (number  of  subscans),  varying  from  ~  0  ’  1 5 
in  the  red  to  ~  0°5  in  the  blue. 

3.2.  Differential  Phases  of  ft  Lyrae 

The  correction  steps  applied  to  the  data  of  the  program  star 
(ft  Lyrae)  to  obtain  the  phase  in  the  Ha  channel  are  presented 
in  Figure  5.  In  order  to  simplify  the  procedures  and  reduce 
the  number  of  calculations,  we  take  advantage  of  the  fact  that 
the  atmospheric  effects  constitute  a  random  and  additive  phase 
term,  so  instead  of  dealing  with  individual  200  ms  subscans  we 
can  work  with  the  average  of  all  visibilities  in  a  scan.  The  left 
panel  of  Figure  5  shows  the  average  phases  of  nine  scans  of  ft 
Lyrae.  In  the  central  panel,  we  show  the  phases  of  the  individual 
scans  after  the  subtraction  of  the  instrumental  phase,  shown  as 
a  dashed  line  in  the  left  panel.  The  last  step  of  the  process  is 


the  correction  of  the  atmospheric  effects  by  fitting  a  quadratic 
function  of  the  form  of  Equation  (3)  to  the  continuum  channels 
for  each  average  scan  and  subtracting  the  resulting  function  from 
all  the  channels.  These  results  are  presented  in  the  right  panel  of 
Figure  5,  where  we  can  see  that  the  resulting  continuum  phases 
are  close  to  zero,  while  the  Ha  ones  are  between  2°  and  5°. 
This  range  of  values  is  due  to  changes  in  the  projected  baseline 
length  over  the  night. 

In  the  procedures  described  above,  the  continuum  phases  of 
ft  Lyrae  were  assumed  to  be  zero.  Nevertheless,  one  should 
expect  to  see  phase  variations  across  the  spectrum,  ft  Lyrae 
is  a  binary  system  with  a  maximum  separation  of  ~  1  mas 
and  Am  ~  0.9  mag,  with  almost  no  color  difference  between 
the  two  components  (Wilson  1974;  Harmanec  2002).  Using 
Equation  (8)  of  Mourard  et  al.  (1992),  we  calculate  that 
observing  this  system  at  maximum  separation  with  a  50  m 
baseline  aligned  with  the  separation  vector  should  show  a 
continuum  phase  0cnt  of  ~  14°  at  the  Ha  channel,  with  a  slope 
of  1°  per  channel.  This  case  corresponds  to  the  longest  baseline 
in  our  data,  but  most  of  our  baselines  were  34  m  long  or  less. 
For  a  30  m  baseline,  </>cnt  at  Ha  is  1°6  with  a  slope  of  0°2  per 
channel. 

Although  the  continuum  phase  component  varies  smoothly 
with  wavenumber  and  can  be  modeled  by  our  atmospheric  phase 
subtraction  method,  arbitrarily  setting  this  component  to  zero 
can,  in  principle,  influence  our  results.  The  complex  visibility 
of  the  Ha  channel  results  from  the  sum  of  the  complex  visibility 
of  the  continuum  and  Ha  emission  line.  Setting  the  continuum 
phase  to  zero  will  change  the  phase  of  the  emission  line  visibility 
by  0Cnt-  For  baselines  of  30  m  or  smaller  this  is  not  an  important 
effect,  given  the  small  continuum  phases.  However,  as  shown 
above,  in  the  case  of  our  longest  baseline  this  can  change  the 
Ha  emission  phase  by  as  much  as  ~  14°.  We  discuss  in  Section 
4  how  we  tested  the  effect  of  this  assumption  on  the  imaging  of 
this  source. 

Figure  6  shows  the  differential  phases  of  ft  Lyrae  on  all 
baselines  observed  on  2005  May  19.  Here  we  see  results  similar 
to  those  from  Figure  5,  with  the  continuum  phases  around 
zero,  and  the  Ha  and  He  i  channels  showing  a  differential 
phase.  The  strongest  Ha  phase  is  found  on  baselines  with 
the  longest  extent  along  the  E-W  direction,  with  very  little 
differential  phase  signals  on  baselines  AN0-AC0  and  AE0- 
AN0.  Comparing  these  results  with  those  from  Mourard  et  al. 
(1992)  and  Harmanec  et  al.  (1996)  we  get  a  consistent  picture 
for  this  system.  Based  on  GI2T  observations,  Harmanec  et  al. 
(1996)  concluded  that  the  orbit  of  the  binary  component  should 
be  roughly  oriented  along  the  east  west  (E-W)  direction,  since 


Figure  5.  Different  correction  steps  applied  to  the  data  of  /3  Lyrae.  The  first  panel  (left)  shows  the  average  phases  of  each  scan  (solid  lines)  for  baseline  E06-AC0, 
obtained  on  2005  May  19.  The  dashed  line  shows  the  instrumental  phases,  obtained  from  the  calibrator  stars,  and  the  dotted  vertical  line  shows  the  location  of  the 
Ho'  channel.  The  middle  panel  shows  the  phases  obtained  after  subtracting  the  instrumental  contribution.  The  right  panel  shows  the  differential  phases  obtained  after 
subtracting  a  quadratic  function,  which  was  fitted  to  the  continuum  points  only  for  each  scan  presented  in  the  middle  panel.  The  uncertainty  in  the  phases  is  of  the 
order  of  0°3  in  the  red,  0°5  at  Ho;,  and  grows  to  1°  in  the  blue. 
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Figure  6.  Residual  phases  of  (3  Lyrae  after  collecting  for  the  instrumental,  atmospheric  and  vacuum  differential  delay  effects.  Each  panel  shows  a  different  baseline, 
observed  on  2005  May  19.  The  individual  lines  represent  individual  scan  averages.  The  vertical  dotted  line  indicates  the  wavenumber  of  Ho;,  while  the  dashed  lines 
show  the  wavenumbers  of  the  He  i  706.5  nm  and  587.6  nm  emission  lines. 


their  observations  used  an  north  south  (N-S)  baseline  and 
did  not  resolve  the  two  stars.  On  the  other  hand,  their  line 
observations  detected  smaller  V2s  for  Ha  and  He  i  indicating 
that  the  line-emitting  region  is  more  extended  along  the  N-S 
direction.  This  result  is  consistent  with  the  spectropolarimetric 
observations  from  (Hoffman  et  al.  1998),  who  found  that  Ha, 
He  i  587.6  nm,  and  ultraviolet  ( X  <  360  nm)  emission  are 
polarized  along  p.a.  ~  74°.  This  result  indicates  that  this 
emission  is  polarized  by  a  structure  perpendicular  to  this 
direction.  Another  result  that  confirms  this  geometry  is  a 
radio  nebula  along  p.a.  156°  ±  4°  detected  by  Umana  et  al. 
(2000). 

Comparing  our  results  to  those  from  Harmanec  et  al.  (1996) 
suggests  a  possible  contradiction,  since  we  do  not  see  a  strong 
differential  phase  on  the  baselines  oriented  closer  to  the  N-S 
direction,  at  least  when  comparing  their  signal  to  that  measured 
on  baselines  oriented  closer  to  the  E-W  direction.  However, 
our  observations  can  be  explained  in  a  way  that  agrees  with 
the  previous  results.  First,  the  weaker  phase  signal  along  the 
N-S  direction  can  be  explained  in  part  as  an  effect  of  the 
baseline  lengths  of  the  two  observations.  Harmanec  et  al.  (1996) 
used  a  baseline  of  51  m,  while  for  our  AE0-AN0  baseline 
we  have  a  longest  projected  length  of  only  ~  25  m  along 
the  N-S  direction.  Our  maximum  baseline  length  along  the 
E-W  direction  is  roughly  twice  that,  which  can  explain  the 
stronger  differential  phase  signal  along  this  direction.  Second, 
we  interpret  the  differential  phases  as  a  displacement  between 


the  Ha  line-emitting  region  and  the  continuum  photocenters. 
Since  the  binary  system  is  expected  to  be  oriented  along  the 
E-W  direction,  this  would  explain  the  stronger  Ha  differential 
phases  on  baselines  aligned  closer  to  this  direction. 

Further  support  for  this  interpretation  is  presented  in  Figure  7, 
where  we  show  the  differential  phases  of  baseline  AW0-E06  on 
2005  May  19  and  2005  May  26.  These  observations  correspond 
to  orbital  phases  of  0.24  and  0.78,  respectively.  Here,  we  can 
see  that  the  Ha  phase  signal  changes  from  negative  to  positive 
between  the  two  dates,  which  can  be  understood  as  the  center  of 
line  emission  and  continuum  photocenter  changing  orientation 
relative  to  each  other. 

We  would  also  like  to  point  out  another  interesting  result 
from  Figure  6,  the  variation  of  the  He  i  differential  phases  near 
k  ~  1.4  and  1.7/rm_1,  in  particular  on  baselines  AE0-AC0 
and  AN0-AC0.  One  can  see  in  this  figure  that  the  He  i  channel 
phases  vary  significantly  from  scan  to  scan,  with  the  sign  of 
the  phase  changing  from  positive  to  negative  in  a  few  scans. 
When  we  compare  this  behavior  to  that  of  the  Ha  phases,  we 
see  that  Ha  does  not  change  sign  on  the  same  night.  This  result 
indicates  that  the  He  i  emitting  region  has  a  structure  different 
from  that  of  the  Ha  region.  This  conclusion  is  supported  by 
the  spectropolarimetric  results  from  Hoffman  et  al.  (1998),  who 
found  that  He  i  587.6  nm  and  He  i  706.5  nm  have  polarization 
properties  different  from  those  of  Ha  and  Hei  667.8  nm. 
This  difference  is  due  to  the  way  these  lines  are  excited,  with 
the  587.6  (k  =  1.70  /nm-1)  and  706.5  nm  ( k  =  1.42 /rm_1) 
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Figure  7.  Comparison  of  the  differential  phases  of  the  longest  baseline  AW0-E06  (53  m)  on  the  nights  of  2005  May  19  (left)  and  2005  May  26  (right).  The  two  nights 
are  separated  by  approximately  half  an  orbital  period. 


being  produced  in  higher  density  regions,  possibly  in  the  inner 
accretion  disk,  consistent  with  the  observed  differences  in  the 
phases  of  these  emission  lines. 

4.  IMAGING 

The  ultimate  goal  of  our  differential  phase  experiment  is 
the  recovery  of  complex  visibilities  of  this  star  in  order  to 
image  the  Ha  and  continuum  emission.  In  order  to  do  this 
we  start  by  combining  the  final  phase  information  described 
above  with  the  calibrated  l/2s  from  the  incoherent  integration 
method  to  calculate  the  corresponding  complex  visibilities  for 
both  continuum  and  line  channels.  Since  these  visibilities  have 
been  corrected  for  instrumental  and  atmospheric  effects,  the 
continuum  channels  can  be  exported  into  a  FITS  file  and  imaged 
using  the  usual  radio  interferometry  techniques.  However,  in 
order  to  image  the  Ha  emission  one  needs  to  take  into  account 
the  fact  that  most  of  the  light  in  this  channel  originates  from  the 
continuum. 

The  correction  for  the  continuum  contribution  to  the  Ha 
channel  is  calculated  using  the  ccnt  value  obtained  from  the 
spectroscopic  measurements.  Rewriting  Equation  (1)  of  Tycner 
et  al.  (2006)  in  a  more  general  form,  we  have  the  following 
relation  for  the  continuum-subtracted  Ha  complex  visibility: 

Vh<x  —  (Fobs  C'cnl  Vent)/ ( 1  Cent)-  (4) 

V0bs  is  the  observed  complex  visibility  in  the  Ha  channel, 
including  the  line  and  continuum  components,  Vcnt  and  VHr/  are 
the  complex  visibilities  of  the  continuum  and  Ha  components, 
and  ccnt  =  0.87  is  the  fractional  contribution  from  the  stellar 
continuum  to  the  total  flux  in  the  Ha  channel  (Section  2.2).  We 
use  the  averaged  visibility  of  the  two  channels  adjacent  to  Ha 
to  estimate  Vcnt  at  the  Ha  channel. 

The  resulting  visibilities  are  converted  into  FITS  files  (one  per 
night)  and  imaged  using  AIPS11  (van  Moorsel  et  al.  1996).  The 
resulting  Ha  images  are  presented  in  Figure  8.  These  images 
were  created  using  pixel  sizes  of  0.5  mas  and  a  hybrid  weighting 
scheme.  We  used  a  combination  between  natural  and  uniform 
weights,  with  robustness  index  0,  which  is  a  good  compromise 
between  resolution  and  sensitivity.  Just  for  reference  we  show 
in  each  panel,  as  a  white  cross,  the  position  of  the  continuum 
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photocenter,  measured  on  continuum  images  restored  from  data 
observed  on  the  same  night.  The  position  of  the  continuum  was 
always  in  the  center  of  the  image,  consistent  with  the  fact  that 
the  continuum  phases  were  set  to  ~  0°  in  our  procedures.  Note 
that  in  most  cases  the  Ha  emission  has  a  shape  similar  to  that 
of  the  restoring  beam,  which  indicates  that  the  emission  is  just 
partially  resolved,  with  a  size  comparable  to  that  of  the  beam. 
We  do  not  find  evidence  of  a  jet  in  these  images,  but  this  does 
not  rule  out  the  presence  of  one,  since  our  lowest  resolution  is 
along  the  N-S  direction,  where  the  jet  is  expected  to  be  found 
(Harmanec  et  al.  1996). 

The  most  important  result  to  be  taken  from  Figure  8  comes 
from  the  comparison  between  the  peaks  of  the  Ha  and  contin¬ 
uum  emission.  We  can  see  in  the  different  panels  that  these  two 
points  do  not  coincide,  and  that  the  displacement  of  the  Ha  pho¬ 
tocenter  changes  relative  to  that  of  the  continuum  as  a  function 
of  orbital  phase.  Since  the  Ha  emission  originates  in  the  disk 
around  the  more  massive  star,  while  the  continuum  photocenter 
is  located  closer  to  the  less  massive  one  (Harmanec  et  al.  1996), 
we  are  effectively  imaging  the  orbit  of  this  system. 

We  also  tested  the  effect  introduced  in  the  Ha  images  of 
assuming  that  the  continuum  phase  on  the  longest  baseline 
is  zero.  To  solve  for  this  effect,  we  added  14°  to  the  phases 
of  both  Ha  and  adjacent  continuum  channels  of  the  longest 
baseline  before  subtracting  the  continuum  contribution  to  the  Ha 
channel.  We  did  not  try  to  apply  this  correction  to  other  baselines 
because  of  the  small  magnitude  of  their  continuum  phases. 
We  followed  the  procedure  described  above  and  obtained  new 
images  using  the  uncorrected  short  baseline  Ha  visibilities  and 
the  corrected  long  baseline  Ha  visibilities.  We  do  not  find  a 
significant  change  relative  to  the  images  obtained  setting  the 
continuum  phase  to  zero.  The  position  of  the  peaks  of  emission 
changed  less  than  0.4er,  indicating  that  setting  the  continuum 
phases  to  zero  was  a  reasonable  approximation  in  this  case. 

5.  THE  ORBIT  OF  THE  Ha  EMISSION  IN  P  FYRAE 

Although  we  barely  resolve  the  stars  in  this  binary,  our  results 
can  be  used  to  derive  some  properties  of  the  ft  Fyrae  system. 
First,  since  this  is  an  eclipsing  binary,  it  is  straightforward  to 
determine  the  orientation  of  the  system  on  the  plane  of  the 
sky.  The  inclination  of  the  orbit  is  86°  (Linnell  et  al.  1998; 
Linnell  2000)  and  its  projection  on  the  sky  can  be  reasonably 
well  approximated  by  a  straight  line.  In  Figure  9,  we  show 
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Figure  8.  Images  of  the  Ha  emission  from  (3  Lyrae,  shown  as  both  grayscale  and  contours.  Each  of  these  figures  corresponds  to  a  different  night,  organized  in  order 
of  increasing  orbital  phase  (top  left:  2005  May  18;  top  right:  2005  May  19;  middle  left:  2005  May  9;  middle  right:  2005  May  24;  bottom  left:  2005  May  26).  The 
cross  shows  the  position  of  the  photocenter  of  the  system,  measured  on  continuum  images.  The  lowest  contour  in  each  figure  corresponds  to  the  3a  noise  level.  The 
reconstructing  beam  is  shown  in  the  bottom  left  corner  of  each  figure. 
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Figure  9.  Positions  of  the  peak  of  Ha  emission  relative  to  the  continuum 
photocenter.  For  each  point  we  show  the  position  uncertainty,  obtained  by  fitting 
a  Gaussian  to  the  image  of  the  source.  Beside  each  point  we  also  indicate  the 
orbital  phase  and  the  date  of  the  observation.  The  dotted  line  was  obtained  by 
fitting  a  least-squares  line  to  the  data.  Since  (3  Lyrae  is  an  eclipsing  binary,  this 
line  is  a  good  approximation  for  the  orientation  of  the  orbit  in  the  plane  of  the 
sky:  Qh«  =  249°. 
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the  position  of  the  Ha  emission  relative  to  the  continuum 
photocenter  for  the  five  nights  of  data  used  in  this  paper.  These 
positions  were  obtained  by  fitting  Gaussians  to  the  images,  and 
the  values  are  given  in  Table  2. 

Figure  9  shows  an  initial  estimate  of  the  position  angle 
of  the  line  of  nodes,  £2Hfy  =  248°8  ±  1°7,  from  fitting  a 
straight  line  to  these  positions.  Note  that  the  point  from  2005 
May  18  has  a  much  higher  uncertainty  in  the  N-S  direction, 
caused  by  the  absence  of  the  northernmost  station,  AN0, 
and  the  resulting  loss  of  resolution  in  this  direction.  We 
fit  an  orbit  for  the  positions  of  Ha  emission  region  around 
the  photocenter  using  instrumental  weighting  (1/a2),  holding 
the  period,  inclination,  and  eccentricity  fixed  and  allowing  the 
semimajor  axis  the  line  of  nodes  £2^,  and  the  time  of 
periastron  passage  just  preceding  our  observations,  7’Ha ,  to  vary. 
We  find  cina  =  0.46  ±  0.03  mas,  £2H„  =  249°  ±  4°,  and 
Tfic  =  ID2453493.77  ±  0.23.  Repeating  this  fit  without  the 
data  point  for  2005  May  9,  when  the  Ha  disk  is  partially  hidden 
by  star  2,  does  not  change  the  results  significantly  (less  than  the 
half  the  quoted  uncertainties). 

The  reliability  of  £2H„  we  obtain  for  the  /l  Lyrae  system  can 
be  checked  by  comparing  with  other  measurements  available  in 
the  literature.  The  radio  observations  presented  by  Umana  et  al. 
(2000)  detected  a  jetlike  nebular  structure  oriented  along  p.a.  = 
156°5  ±  4°.  This  nebula  is  perpendicular,  within  uncertainties, 
to  the  orbit  derived  here,  agreeing  with  the  expectation  that 
jets  are  launched  perpendicular  to  their  disks.  Further  evidence 
confirming  the  orbital  orientation  of  this  system  comes  from  the 
comparison  of  our  results  with  the  spectropolarimetric  results 
of  Hoffman  et  al.  (1998).  They  predict  that  the  jets  should  be 
oriented  along  p.a.  ~  163°5,  similar  to  that  found  by  Umana 
et  al.  (2000)  and  roughly  perpendicular  to  the  orbital  orientation 
we  found. 

Finally,  we  compare  our  measurements  of  the  displacement 
between  the  Ha  and  continuum  photocenters  to  predictions 
based  on  parameters  of  the  system  derived  by  others  (e.g., 
Harmanec  et  al.  1996;  Harmanec  2002;  Linnell  2000).  We  use 
the  reference  system  of  Harmanec  et  al.  (1996),  where  the  orbit 
is  in  the  X  -  Y  plane,  with  the  origin  at  the  most  massive 
star  (the  one  that  is  currently  gaining  mass)  and  the  X-axis 


Figure  10.  Ha  emission  positions  compared  to  model  orbits.  The  solid  line 
shows  the  model  orbit  based  on  the  parameters  discussed  in  Section  5.  The 
dotted  line  shows  the  same  orbit,  but  with  an  assumed  distance  of  309  pc,  la 
greater  than  the  Hipparcos  distance.  The  lines  connecting  the  emission  positions 
to  the  orbit  indicate  the  model  position  along  the  orbit  for  each  night. 


defined  by  the  line  joining  the  two  stars.  The  analysis  of  the 
radial  velocities  of  this  system  shows  that  the  eccentricity  of  the 
orbit  is  0  (Harmanec  &  Scholz  1993).  According  to  Harmanec 
(2002),  the  separation  between  the  centers  of  the  two  stars  is 
58.52  Rq.  For  A m  =  0.9  mag  the  position  of  the  photocenter 
is  at  X  =  40.64  RQ  and  Y  =  0.  Based  on  the  analysis  of  the 
velocity  curve  of  the  Ha  emission,  Harmanec  et  al.  (1996)  found 
that  this  emission  line  originates  in  the  disk  at  X  —  4.12  Rq, 
Y  =  4.97  Rq,  a  region  offset  from  the  center  of  the  system, 
possibly  related  to  the  point  of  impact  of  the  gas  in  the  disk. 

We  used  these  parameters  to  calculate  the  model  orbits  of 
the  Ha  and  continuum  photocenters,  and  converted  them  to 
positions  in  the  sky  using  the  inclination  derived  by  Linnell  et  al. 
(1998;  i  —  86°),  the  orientation  derived  by  us  (£l[\„  =  249°), 
and  a  distance  of  270  pc.  The  resulting  Ha  orbit  relative  to  the 
continuum  photocenter  is  shown  in  Figure  10  as  a  solid  line. 

The  Ha  orbit,  at  «Hc,  =  0.46  ±  0.01  mas,  is  smaller  than  the 
model  Ha  orbit  (flmodei  =  0.63  mas)  by  26%.  There  are  several 
possible  explanations.  First,  the  system  may  be  more  distant 
than  270  pc.  Increasing  the  distance  by  39  pc  (la)  would  reduce 
the  model  orbit  size  by  14%.  This  suggestion  is  supported  by 
a  recent  new  reduction  of  the  Hipparcos  data  (van  Leeuwen 
2007),  which  finds  a  distance  of  315+ll/  pc  for  fi  Lyrae.  Second, 
the  magnitude  difference  may  be  smaller  than  0.9  mag;  if  Am  — 
0.8  mag,  the  photocenter  is  closer  to  star  1,  hence  closer  to  the 
Ha  emission,  and  the  Ha  orbit  is  reduced  by  ~  3%.  However, 
we  will  show  in  the  next  section  that  the  magnitude  difference 
between  the  stars  is  likely  to  be  higher  than  0.9  mag. 

Finally,  an  important  source  of  error  in  the  size  of  the  mea¬ 
sured  Ha  orbit  is  the  correction  for  the  continuum  contribution 
to  the  Ha  channel.  Factors  that  contribute  to  errors  in  this  cor¬ 
rection  are  the  uncertainty  in  the  width  of  the  channel  and  in 
the  EW  of  the  Ha  absorption  component.  In  this  paper,  we  used 
the  Cote  &  Waters  (1987)  calibration  to  correct  for  the  Ha  ab¬ 
sorption  in  the  stellar  spectrum  to  obtain  ccnt  =  0.87  ±  0.01  as 
the  fractional  contribution  from  the  stellar  photosphere  to  the 
flux  in  the  Ha  channel.  However,  their  calibration  is  based  on 
main-sequence  stars  and  may  overestimate  the  correction  for 
(1  Lyrae,  where  the  brightest  component  is  a  B6-8II  star.  We 
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Table  2 

Displacement  of  the  H a  and  Continuum  Photocenters 


UT  Date  Orbital  Phase  ccnt  =  0.87 


Cent  =  0.91 


R.A. (Ha -Continuum)  Decl. (Ha -Continuum)  R.A. (Ha -Continuum)  Decl. (Ha -Continuum) 

(mas)  (mas)  (mas)  (mas) 


2005  May  9 

0.47 

-0.13 

± 

0.03 

-0.02 

± 

0.05 

-0.33 

± 

0.04 

-0.13 

± 

0.07 

2005  May  18 

0.17 

-0.44 

± 

0.05 

-0.05 

it 

0.16 

-0.75 

± 

0.04 

-0.06 

it 

0.20 

2005  May  19 

0.24 

-0.44 

± 

0.03 

-0.16 

± 

0.04 

-0.76 

± 

0.03 

-0.25 

± 

0.05 

2005  May  24 

0.63 

0.29 

± 

0.03 

0.15 

± 

0.05 

0.52 

± 

0.04 

0.32 

± 

0.07 

2005  May  26 

0.78 

0.39 

it 

0.02 

0.18 

± 

0.03 

0.73 

± 

0.04 

0.33 

± 

0.06 

Notes.  Columns  3  and  4  show  the  measurements  obtained  from  images  created  using  ccnt  =  0.87,  while  Columns  5  and  6  show  the 
measurements  from  images  created  using  ccnt  =  0.91.  The  orbital  phases  were  calculated  based  on  primary  eclipse  ephemerides 
from  Kreiner  (2004). 


Figure  11.  Ha  emission  positions  derived  using  a  larger  correction  for  the 
continuum  contribution  to  the  Ha  channel  (cCOnt  =  0.91).  The  solid  line  shows 
the  model  orbit  for  a  distance  of  270  pc. 
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Figure  12.  Separation  of  the  two  binary  components,  obtained  by  fitting  the 
continuum  V2  data.  For  each  point  we  show  the  uncertainty  in  the  position  of 
the  separation  of  the  two  sources.  Beside  each  point  we  also  indicate  the  orbital 
phase  and  the  date  of  the  observation.  The  dotted  line  was  obtained  by  fitting 
a  least  squares  line  to  the  data,  and  approximates  the  orientation  of  the  orbit  in 
the  plane  of  the  sky  (£lbin  =  256°). 
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tested  the  dependence  of  the  orbit  size  on  this  correction  by 
using  the  extreme  assumption  that  the  EW  correction  for  the 
underlying  absorption  to  Ha  is  0,  which  gives  ccnt  =  0.91.  Us¬ 
ing  this  value,  we  obtained  new  images  and  remeasured  the  Ha 
emission  positions  with  respect  to  the  continuum  photocenter 
(Table  2).  The  result  was  to  increase  the  separations  by  a  factor 
of  ~  1.8,  which  increases  a to  ~  0.85  mas,  or  ~  30%  larger 
than  amodei  (Figure  11).  We  can  infer  that  even  alcr  increase  in 
Cent  would  increase  aua  to  ~  0.545  mas. 

6.  COMPARISON  BETWEEN  DIFFERENTIAL  PHASE 
AND  V2  RESULTS 

As  a  way  to  check  the  Ha  orbit  obtained  with  the  differential 
phase  technique,  we  use  the  calibrated  continuum  V2  and 
closure  phases  to  determine  the  separation  p,  position  angle 
0,  and  magnitude  difference  of  the  binary  for  each  night.  The 
Harmanec  et  al.  (1996)  model  of  the  locations  of  the  stellar 
components  and  the  Ha  emission  within  the  /3  Lyrae  system 
then  allow  us  to  determine  whether  the  binary  and  the  Ha  orbits 
are  consistent  with  one  another. 

To  determine  the  binary  orbit,  we  fit  p,  6,  AV,  and  A R  to 
each  night’s  V2  and  closure  phase  data.  We  did  a  grid  search  for 
the  parameters  that  minimize  /  2 ■  The  range  of  parameters  used 
was  0.5  <  p  <  2.0  mas  in  steps  of  0.01  mas,  0°  <  9  <  360°  in 


steps  of  CL01,  0.7  <  AV  <1.9  mag,  and  0.7  <  A R  <1.9  mag 
in  steps  of  0.01  mag.  We  obtained  the  values  given  in  Table  3 
and  plotted  in  Figure  12.  Note  that  we  could  not  converge  on 
a  fit  for  the  night  2005  May  09,  when  a  significant  part  of  star 
1  is  hidden  behind  the  donor  star.  The  best-fit  orbit  to  these 
points  has  «bin  =  0.99  ±  0.03  mas,  fibin  =  255°  ±  7°,  and 
7],in  =  ID2453493.47  ±  0.34,  where  the  subscript  “bin”  refers 
to  the  binary  orbit,  as  distinct  from  the  Ha  orbit.  The  orbital 
parameters  obtained  here  as  well  as  those  obtained  fitting  the 
Ha  measurements  are  given  in  Table  4. 

The  value  for  ahm  is  very  close  to  the  predicted  value  of  1.01 
mas.  The  uncertainty  in  the  orbital  positions  is  large,  in  part 
because  the  baselines  we  used  do  not  fully  resolve  this  system, 
as  can  be  seen  in  Figure  13.  Here  we  present,  for  each  night,  the 
V 2  of  one  scan  taken  with  the  longest  baseline  (AW0-E06, 53  m) 
within  30  minutes  of  meridian  crossing,  when  the  resolution  for 
this  baseline  is  highest.  This  figure  shows  some  of  the  limitations 
of  our  measurements.  As  one  would  expect,  we  observe  higher 
V 2  on  the  night  of  May  24,  when  the  system  is  close  to  eclipse, 
and  lower  V2  on  May  19  and  May  26,  when  the  separation 
of  the  stars  is  close  to  its  maximum.  However,  our  baselines 
are  not  long  enough  to  show  the  full  characteristic  cosine  wave 
signature  of  a  binary,  which  increases  the  uncertainty  in  the 
determination  of  the  separation  and  magnitude  difference  of  the 
two  stars.  Future  work,  using  longer  baselines,  will  be  able  to 
reduce  the  uncertainties  in  these  values. 
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Table  3 

Binary  Separations 


UT 

Date 

Orbital 

Phase 

P 

(mas) 

e 

(deg) 

AMaj 

(mas) 

^Min 

(mas) 

<t> 

(deg) 

AV 

(mag) 

A  R 
(mag) 

2005  May  18 

0.17 

0.93 

263.1 

1.09 

0.19 

177.9 

1.38  ±0.1 

1.20  ±0.1 

2005  May  19 

0.24 

0.95 

266.1 

0.52 

0.21 

2.7 

1.23  ±  0.1 

1.20  ±0.1 

2005  May  24 

0.63 

0.85 

74.4 

0.51 

0.23 

3.8 

1.24  ±0.1 

1.20  ±0.1 

2005  May  26 

0.78 

0.94 

87.7 

0.52 

0.21 

1.9 

1.34  ±0.1 

1.18  ±  0.1 

Notes.  Column  1 :  date  of  the  observations;  Column  2:  orbital  phase;  Columns  3  and  4:  binary  separation  and 
position  angle;  Columns  5-7:  position  uncertainty  ellipse  major  axis,  minor  axis,  and  position  angle;  Columns  8 
and  9:  magnitude  difference  between  the  two  sources. 


Table  4 

Orbital  Elements 


Data 

Value 

aHa  (mas)a 

0.46  ±  0.03 

aKa  (mas)b 

0.80  ±  0.03 

t/bin  (mas) 

0.99  ±  0.03 

i (deg)c 

86 

£2Ha  (deg) 

249  ±4 

£2bi„  (deg) 

255  ±  7 

0 

o)  (deg) 

270 

7h„  UD) 

2453493.77  ±  0.23 

7bi„  (JD) 

2453493.47  ±  0.34 

P  (days)e 

12.9421 

Notes.  rHc  and  7bin  correspond  to  the  fitted  time  of 
periastron  passage  just  preceding  our  observations. 
a  Value  obtained  from  images  created  assuming  ccnt  = 

0.87. 

b  Value  obtained  from  images  created  assuming  ccnt  = 

0.91. 

c  Data  from  Linnell  et  al.  (1998). 
cl  Data  from  Harmanec  &  Scholz  (1993). 
e  Data  from  Kreiner  (2004). 

To  compare  the  Her  orbit  with  the  binary  orbit,  we  start  with 
the  model  discussed  in  Section  5:  the  two  stars  are  separated 
by  58.52  Re\  the  continuum  photocenter  is  located  43.96  Rq 
from  the  more  massive  star  (star  1),  calculated  using  the  average 
A R  =  1.20  mag  from  Table  3;  and  the  Her  emission  is  39.84  Rq 
from  the  photocenter,  close  to  star  1  but  4.97  RQ  off  the  star- 
to-star  line  (Harmanec  et  al.  1996).  This  model  implies  that 
the  Her  orbit  around  the  continuum  photocenter  should  have  a 
semimajor  axis  68%  as  big  as  that  of  the  binary  orbit,  and  it 
should  lag  star  1  in  the  binary  orbit  by  ~  0.02  in  orbital  phase, 
or  0.25  days. 

In  Section  5,  we  found  that  aH»  =  0.46  ±  0.03  mas,  con¬ 
siderably  smaller  than  68%  of  the  observed  c/hm  =  0.99  ± 
0.12  mas.  However,  as  also  discussed  in  Section  5,  the  formal 
0.01  uncertainty  in  ccnt  corresponds  to  a  ~  0.085  mas  uncer¬ 
tainty  in  aHo,  which  is  larger  still  if  we  take  into  account  the 
further  uncertainty  of  applying  the  Cote  &  Waters  (1987)  cali¬ 
bration  to  the  luminosity  class  II  primary  of  ft  Lyrae.  Given  the 
uncertainties,  aHol  and  ahln  are  not  inconsistent. 

The  orientations  of  the  orbits  derived  from  the  two  methods 
are  quite  consistent.  Our  orbital  fits  give  £i[\u  =  249°  ±  4°  and 
£2bin  =  255°  ±7°,  which  agree  well  within  the  uncertainties.  The 
phase  lag  7j \(1  —  7bjn  between  our  calculated  time  of  periastron 
passage  for  the  Ha  and  binary  orbits,  0.30  ±  0.41  days  or 
0.023  ±  0.032  in  orbital  phase,  is  also  in  very  good  agreement 
with  the  model  predictions. 

Another  important  result  obtained  from  the  analysis  presented 
here  is  the  magnitude  differences  between  the  two  stars.  The 
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Figure  13.  Calibrated  V2  data  from  J3  Lyrae  observed  with  the  longest  baseline, 
AW0-E06  (53  m).  The  respective  nights  and  orbital  phases  are  indicated  inside 
the  panels.  In  order  to  make  the  visualization  of  the  plots  easier  each  panel 
shows  only  one  scan,  corresponding  to  the  one  observed  within  30  minutes 
of  local  meridian  crossing.  The  solid  lines  show  the  best  fitting  model  for  the 
corresponding  night. 

value  used  in  Section  5  (Am  =  0.9  mag)  was  derived  by 
Wilson  (1974)  using  simple  models.  Based  on  a  more  detailed 
model,  which  fits  the  light  curve  at  several  bands  simultaneously, 
Linnell  et  al.  (1998)  found  AT  =  1.8  mag  at  orbital  phase  0.25. 
Our  results  show  significantly  smaller  values,  AV  —  1.30  ± 
0.1  mag  and  A R  =  1.20  ±0.1  mag,  which  will  have  important 
implications  for  the  models  of  this  binary  system. 

As  a  last  check  we  used  the  method  described  by  Tycner  et  al. 
(2005,  2006)  to  determine  the  size  of  the  Ha  disk.  We  used 
the  calibrated  Ha  V2,  corrected  to  eliminate  the  continuum 
contribution  from  the  binary  star  to  this  channel.  Following  the 
results  from  Tycner  et  al.  (2006)  we  decided  to  fit  the  Ha  V2  with 
a  Gaussian  model,  which  results  in  a  disk  with  HWHM  =  0.60  ± 
0. 10  mas  (Figure  14).  The  Gaussian  half  width  at  half-maximum 
(HWHM)  is  comparable  to  the  Roche  lobe  radius  of  ~  0.52  mas 
(30.3  Rq  at  270  pc;  Harmanec  2002),  and  consistent  with  the 
orbit  obtained  using  the  two  methods  described  above. 

7.  SUMMARY 

In  this  paper,  we  presented  the  development  of  a  differential 
phase  technique  using  NPOI  data.  We  described  the  methods 
used  to  correct  the  instrumental  and  atmospheric  effects  on 
the  phases  of  ft  Lyrae,  which  allowed  us  to  recover  the 
complex  visibility  of  the  Ha  channel.  We  also  described 
the  procedure  used  to  correct  the  continuum  contribution  to 
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Figure  14.  Calibrated  V2  data  for  the  Hey  channel  as  a  function  of  radial  spatial  frequency,  ( u 2  +  tr)1/2.  The  dashed  line  shows  the  continuum  V2  after  the  removal  of 
the  binary  signature.  The  dotted  line  shows  the  best  elliptical  Gaussian  model  fitted  to  the  data. 


the  Ha  channel  and  successfully  imaged  the  line  emission  of 
this  system  using  standard  radio  interferometric  techniques.  We 
found  that  the  major  axis  is  oriented  along  £2  =  249°,  consistent 
with  previous  spectropolarimetric  measurements  and,  radio  and 
optical  interferometry  results.  From  the  comparison  between  the 
observed  position  of  Ha  relative  to  the  continuum  photocenter 
and  values  obtained  from  models,  we  find  that  our  measurements 
indicate  a  semimajor  axis  for  the  Ha  orbit  that  is  smaller  than 
that  predicted  by  the  models.  We  suggest  that  the  most  likely 
causes  for  this  discrepancy  are  errors  in  the  distance  and  in  the 
correction  for  the  continuum  contribution  to  the  Ha  channel.  We 
found  that  the  results,  obtained  using  the  technique  developed 
here,  are  consistent  with  those  obtained  from  a  more  traditional 
analysis  that  uses  only  V 2  and  closure  phase  measurements. 
Also,  based  on  the  V 2  analysis,  we  were  able  to  determine 
the  size  of  the  Ha  disk  and  the  magnitude  difference  between 
the  two  stars,  which  is  an  important  input  to  models  of  this 
system. 

The  technique  presented  here  represents  a  major  step  in 
dealing  with  the  problem  of  the  loss  of  phase  information 
in  optical  interferometric  observations.  At  visual  wavelengths, 
the  phase  fluctuations  caused  by  atmospheric  turbulence  are 
both  large  and  rapid,  badly  corrupting  the  measured  fringe 
phases.  Consequently,  most  optical  interferometric  imaging  has 
been  carried  out  using  the  square  of  the  fringe  amplitude, 
V2,  and  the  closure  phase,  in  which  the  atmospheric  phase 
fluctuations  cancel.  By  making  some  assumptions  about  the 
structure  of  the  source,  our  technique  allows  us  to  recover 
the  fringe  phase  information,  and  thus  to  produce  complex 
visibilities.  With  these  complex  visibilities  in  hand,  we  can 


apply  the  image  reconstruction  techniques  that  are  employed 
in  radio  interferometry. 

In  the  case  of  /3  Lyrae,  this  technique  allowed  us  to  image  the 
Ha  emission  and  determine  its  location  and  motion  with  respect 
to  the  photocenter  of  the  system,  and  thus  to  determine  its  orbit. 
We  also  found  that  the  phases  related  to  the  He  i  emission  have  a 
behavior  different  from  that  of  the  Ha,  suggesting  that  they  are 
formed  in  different  regions.  As  discussed  in  the  introduction, 
there  is  evidence  that  the  Ha  emission  of  this  system  may  origi¬ 
nate  in  a  jet  (Harmanec  et  al.  1996;  Hoffman  et  al.  1998;  Umana 
et  al.  2000).  We  do  not  see  evidence  for  a  jet  in  our  images. 
However,  this  may  be  due  to  the  lower  resolution  along  the  N-S 
direction,  which  is  the  direction  along  which  the  jet  is  expected 
to  be  oriented.  Future  observations  with  longer  baselines  will 
allow  us  to  put  better  constrains  on  the  existence  and  size  of  a  jet. 

Besides  p>  Lyrae  we  are  also  applying  this  technique  to 
other  Be  systems.  This  technique  will  allow  us  to  image  their 
circumstellar  disks  and  map  their  structures,  and  study  in  better 
detail  the  accretion  flow  and  effects  such  as  the  ionization  of  the 
disk  by  a  binary  companion. 

The  work  done  with  the  NPOI  was  performed  through  a 
collaboration  between  the  Naval  Research  Lab  and  the  US  Naval 
Observatory,  in  association  with  Lowell  Observatory,  and  was 
funded  by  the  Office  of  Naval  Research  and  the  Oceanographer 
of  the  Navy.  We  thank  the  NPOI  staff  for  the  careful  observations 
that  contributed  to  this  work.  C.T.  thanks  Lowell  Observatory  for 
the  generous  time  allocation  on  the  John  S.  Hall  1.1m  telescope 
and  thanks  Wes  Lockwood  and  Jeffrey  Hall  for  supporting  the 
Be  star  project  on  the  Solar-Stellar  Spectrograph.  This  research 
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